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Abstract. We present measurements of growth rates of the (0001) facet sm-face of ice as a 
function of water vapor supersaturation over the temperature range —2 > T > —40 C. From these 
data we infer the temperature dependence of premelting on the basal surface and the effects of 
premelting on the ice growth dynamics. Over this entire temperature range the growth was consistent 
with a simple 2D nucleation model, allowing a measurement of the critical supersaturation cfq(T) as 
a function of temperature. We find that the 2D nucleation barrier is substantially diminished when 
the premelted layer is partially developed, as indicated by a reduced ctq, while the barrier is higher 
both when the premelted layer is fully absent or fully developed. 

1 Introduction 

It has long been long suspected that premelting in ice plays an important role in the growth dy- 
namics of ice crystals from water vapor [U O |3] . Although ice is a monomolecular crystal with a 
simple hexagonal structure under normal atmospheric conditions, ice crystals forming from water 
vapor exhibit an exceedingly rich spectrum of plate-like and columnar morphologies as a function 
of temperature and supersaturation over the temperature range > T > —30 C [3]. Since the pre- 
melted layer in ice develops over this same temperature range [IJ [5j [6] , the prevailing thinking holds 
that temperature-dependent effects of premelting on ice crystal growth are responsible for the ob- 
served morphological complexities, together with instabilities arising from diffusion-limited growth 
and other effects [3]. To date, however, this long- held hypothesis has remained largely unsupported 
by solid experimental evidence. 

We sought to shed light on this problem by making precise measurements of the growth rates of 
small faceted ice crystals from water vapor under carefully controlled conditions, in order to better 
quantify the intrinsic ice growth behavior. To this end we measured growth rates of (0001) ice 
surfaces as a function of water vapor supersaturation a and temperature T over the temperature 
range —2 > T > —40 C. Our measurements were made at low background pressure to reduce the 
effects of particle diffusion, so the growth was mainly limited by attachment kinetics. 

As we will show below, our measurements were all consistent with a simple 2D nucleation model 
of crystal growth, and from our data we extracted the critical supersaturation ctq (T) as a function of 
temperature. The function cro(r) showed an interesting behavior from which we were able to infer 
the onset and development of premelting on the ice surface as well as how premelting affects crystal 
growth in ice. 

2 Ice Crystal Growth Measurements 

The goal of our ice growth experiments was to examine the growth of individual ice crystals in a 
carefully controlled environment, and an idealized schematic diagram of our experimental set-up is 
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Figure 1: An idealized schematic of our experimental set-up. The top surface is an ice reservoir at 
temperature Tm that supplied water vapor for a test crystal resting on a substrate at temperature 
Tsubst- When Tm > Tgubst, growth rates were determined by measuring the size and thickness of 
the test crystal as a function of time. 



shown in Figure [T] The top surface of the experimental chamber was a thermal conductor with a 
uniform temperature Tm, and its inside surface was covered with a layer of ice crystals that made 
up an ice reservoir. At the beginning of each measurement, a single test crystal was placed near 
the center of the bottom substrate surface held at temperature Tsubst- The walls were thermally 
insulating, and the vertical spacing from the top to the bottom of the chamber was 1.0 mm. The 
temperature difference AT — Tjn ~ Ts^b^t determined the effective supersaturation seen by the test 
crystal. During the course of a measurement we increased AT and observed the size and thickness 
of the test crystal as a function of time, and from this we extracted crystal growth velocities as a 
function of supersaturation. A detailed hardware description is provided in [7]. Some important 
aspects of the experiment included: 

1) Our test crystals were small, thin plates, typically < 5 /im thick and < 50 \xm in diameter. 
This was important to reduce the effects of particle and heat diffusion, as described quantitatively 
inC]. 

2) We used crystals with one basal facet in contact with the substrate, measuring the perpendicu- 
lar growth of the opposite (0001) surface using broad-band interferometry as described in 7 . Growth 
of the prism facets was affected by substrate interactions, particularly at low supersaturations, so 
these measurements were discarded. 

3) We used only crystals with simple morphologies and well formed facets, and each crystal was 
discarded after growth. Evaporating and regrowing crystals was found to result in generally lower 
quality data. 

4) Our test crystals were freshly made in a clean environment and transported within minutes 
to our test chamber with minimal processing, as described in [7]. 

With these precautions, and exercising considerable care to create a stable growth environment 
with precisely known supersaturations, we were able to obtain quite satisfactory data, and an ex- 
ample is shown in Figure [51 Here we have plotted the effective attachment coefficient a, which is 
derived from the perpendicular growth velocity v using v = avkin<^, where a is the supersaturation 
far from the crystal and Vkin is a kinetic velocity [3) [7] . 
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Figure 2: Measurements of the growth of the basal facets of two ice crystals at -15 C, shown as 
the effective condensation coefficient a as a function of supersaturation <j far from the crystal. One 
crystal (dots) was grown in a background pressure of air at 20 Torr and the other (squares) was grown 
in a background pressure of 740 Torr. The low-pressure crystal shows mainly kinetics-limited growth, 
while the growth at high pressure is mainly limited by particle diffusion when the supersaturation 
is high. Fit lines are described in the text. 

For our main data set, each crystal was grown at a pressure near 20 Torr, where the growth was 
predominantly limited by attachment kinetics, at least for small crystals at low supersaturations. 
For essentially all our data, we found that the basal growth was well described by a 2D nucleation 
model [9l, and to describe the growth we adopted a simplified parameterization of the intrinsic 
attachment coefhcient aintrinsicio',T) ~ Aexp{—ao/a) where A and (Tq ^^re parameters that may 
depend on temperature but not on supersaturation. 

At high supersaturations our data were distorted by the effects of particle and heat diffusion, so 
we modeled our data by fitting to the functional form 



as described in 7 . To produce more stable fits in our analysis, we adjusted ctq a-nd a fa while 
keeping the scale factor A fixed, and we chose A — 1 to give the physically reasonable result that 
^intrinsic — >■ 1 as (7 bccomcs large. In Figure [H for example, the fit parameters are (A, ijoi ct/it) 
= (1,2.3,0.15) and (1,2.5,0.0075) for the low-pressure and high-pressure crystals, respectively. As 
described in [7] , the value of a fu depended mainly on the background pressure as well as the size and 
thickness of each test crystal, so this parameter was of little physical interest in our later analysis; 



Aexp (-o-q/ct) a fit 



(1) 



Aexp {-(Jo/ a) -I- a/^ 



3 



the main focus in this experiment was on measuring the critical supersaturation ctq as a function of 
temperature. 

Figure |3] shows data from a typical six-hour run. Although each crystal took only about ten 
minutes to measure, finding suitable crystals added time, plus there was additional procedural 
overhead, so the experiment ended up taking about 45 minutes per crystal. For each measurement, 
like those shown in Figure [31 we generated curves bracketing the fit curves as a visual means of 
producing an error estimate Sao for each measured ctq- The dotted lines in Figure |3] show the fit 
curves for each data set except with ctq replaced by (Tq ± Sao in the fitting function Equation [T] In 
total, 102 crystals were measured and fit in this way to produce our primary reduced data set. 

We combined the individual i7o,i measurements at each temperature by weighting each measure- 
ment with our estimated Sao,i for that measurement, thus producing a weighted mean {ao) at each 
temperature along with an uncertainty estimate S (ao) for the mean [8] 

Y.'^o,iSaQ'^ 



S{aof = 



iNeff~l)E3<Jo^^ 



where N^ff — ^'^o'i) /'^'^cTmin the effective number of points in the sample at each temperature. 

We examined our data fitting and analysis practices carefully and reached several conclusions: 
1) Nearly all our data were well fit by Equation [1] consistent with a simple 2D-nucleation model 
distorted by particle diffusion; 2) Roughly ten percent of the crystals sampled grew very rapidly, 
with essentially no nucleation barrier These crystals may have had dislocations, and they were 
discarded from our data. These crystals were obviously different from the norm, however, and we do 
not believe this practice was detrimental to our analysis; 3) Our fits were quite robust to the critical 
supersaturation ao, since this term was determined mainly by the growth behavior at low a, while 
diffusion predominantly affected the high-cr growth. Since our data were taken at low pressure, the 
inclusion of a fa mainly helped fit the higher-tr part of the each curve, which was useful for producing 
a better global fit. 

Averaging our data at each temperature yielded the final measurement of the critical supersatu- 
ration as a function of temperature ao (T) shown in Figure 31 which is the principal result from our 
experiment. Here the current set of measurements includes only data points for T > —20 C; the 
last two data points in Figure HI were taken from [10) . which was a previous version of the present 
experiment. 

Since the crystal growth we observed was everywhere well described by a 2D nucleation model, 
the measured critical supersaturation ao can be used to calculate the step free energy /3 using 

^ 3fc2r2 

where 0,2 is the area of a molecule on the surface, a relation that comes from classical 2D nucleation 
theory [9l|T0]. A plot of /3{T) from our data is shown in Figure |5l We note from the scale on the 
right side of Figure El that /3(T) is much smaller than /3o = 7a « 3.5 x 10^^^ J/nr, the product of 
the surface energy 7 = 0.11 J/m^ of the ice/vapor interface and the nominal molecular step height 
a = 0.32 nm, which is an upper limit on the step energy [lOj . Since the step energy is an equilibrium 
quantity, its calculation using perturbation techniques or molecular dynamics simulations may be a 
tractable problem, so this is an area for further research. 
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Figure 3: Eight separate data sets showing the growth of crystals at -12 C. Individual measurements 
of a (cr) have been offset by multiples of 0.05 for clarity; the amount of each offset can be seen by 
noting that a (0) = 0. Solid lines drawn through each data set show fits to the functional form 
described by Equation [1] and each curve is labeled with the fit critical supersaturation oq. 
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Figure 4; Combined measurements of the critical supersaturation (Tq as a function of temperature, 
as described in the text. 



2.1 The Prefactor A 

A remaining question concerns our choice to use a constant prefactor A in our fitting procedure, and 
our choice of the numerical value of ^ = 1 . We explored this question at some length by examining 
the behavior of our data, and an example is shown in Figure |6l The eight crystals displayed in this 
plot were chosen because they exhibited well-behaved growth profiles a{a) with minimal growth 
perturbations from diffusion, and the residual diffusion effects were fit and corrected before plotting 
the data. If the growth satisfies a nucleation model with amtrinsic = ^exp(— cto/ct), then by plotting 
the data as shown in Figure |6l the extrapolated data should intercept the (t~^ — axis at the point 
where a = A. The lines drawn through the points used the individual CTq fit values as described 
above, which assumed A= 1. 

By examining our data in this fashion, both with and without diffusion corrections, it became 
clear that a prefactor o{ A = 1 was a good approximation for the full range of temperatures we 
measured, at least to an accuracy of about a factor of two. It was difficult to make an accurate 
measurement of A, however, because this involved extrapolating our data, and such a measurement 
was somewhat affected by diffusion effects. Including A as a free parameter in our fitting procedure 
for individual crystals gave the fits too many degrees of freedom and produced unstable results. We 
therefore chose to fix A, and we chose A — 1 from what was essentially a global fit to the data at 
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Figure 5: Step free energy /3(T) as a function of temperature, calculated from the critical supersat- 
uration (To(T') according to classical nucleation theory, as described in the text. The scale on the 
right side of the graph shows (3(T) normalized by /3o, the product of the surface energy and the step 
height [lOj. 

all temperatures. As mentioned above, this value also gives the physically reasonable result that 
aintrinsic -> 1 as fj becomes large. 

To see whether our choice of A had a strong effect on our final results, we performed an exercise 
of changing from A = 1 to A = 0.5 and reanalyzing our data. For this we refit each crystal to 
Equation [1] using A — 0.5, generating a new ctq and afit for each crystal, and from the new data set 
we recombined the data at each temperature, following the procedure outlined above. The results 
are shown in Figure [71 The points, error bars, and the solid line in this figure are identical to Figure 
m showing the original analysis with A = I. The new analysis with A = 0.5 produced a new set 
of points shifted downward, and these are shown as a dashed line in the figure; the error bars are 
essentially unchanged from the original analysis. 

Although A ^ 1 gives a somewhat better global fit to our data, a value oi A = 0.5 is not outside 
the range of potential systematic errors in our experiment. From our analysis including plots like 
that shown in Figure [HI we believe that A = 1 is probably accurate to within a factor of two or 
better. However we see from Figure [7] that changing A by a factor of two produced only a modest 
shift in our final measurement of ao{T). 
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Figure 6: A plot of the measured attachment coefficient a as a function of the inverse supersaturation 
1/a for a representative sample of two crystals at each of four different temperatures, as labeled in 
the plot. Fit lines are drawn to show that the extrapolated data are all consistent with a prefactor 
oi A= 1, as described in the text. 

3 Interpretation and Discussion 

In summary, we measured growth velocities of the (0001) surface of ice as a function of water vapor 
supersaturation over the temperature range —2 < T < —40. From the measured growth velocities 
we extracted the intrinsic attachment coefficient aintrinsic using v = aintrinsicVkin<^ sur f , where here 
o'surf is the supersaturation immediately above the ice surface. Over the entire temperature range, 
and for all supersaturations measured, our data were well described by a simple 2D nucleation 
model using an attachment coefficient aintrinsic = Aexp{—ao/a). A global fit to the data yielded 
the constant prefactor A = I, accurate to better than a factor of two over our entire temperature 
range. Fits to the data at each temperature yielded the measured critical supersaturation (Jq(T) 
shown in Figure 01 

Referring to Figure SI let us now examine several different temperature regions and discuss the 
behavior of uo(T) in detail. 

Region III: [—13 > T > —40 C]. In this temperature region we see a simple power-law behavior 
in cro(r), which was observed previously in [10]. From the data in Region III we have o'o(T) « 
0.009T2 percent and (3{T) k, 1.2 x lO^^^ j/j^, where for both these T is in degrees C. 

We believe premelting is essentially absent in this temperature region, meaning that whatever 
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Figure 7: Combined measurements of the critical supersaturation ctq as a function of temperature, 
along with a reanalysis using the prefactor A — 0.5. The points, error bars, and solid line are the 
same as those shown in Figure S) The dashed line shows how the results shifted downward when 
the analysis was redone using A = 0.5. 

premelting effects do exist are too small to affect the ice growth dynamics. We have no theoretical 
explanation for the observed temperature dependence of a(){T). As the temperature increases, 
however, we might expect that surface restructuring would be more likely to smooth out the edge 
of a 2D island, thus lowering the step energy. Quantitative calculation of this effect is an area for 
further research. 

The peak at T = —12 C. We examined the small peak in (To at this temperature rather carefully, 
and it does not appear to be a statistical fluctuation, nor could we find any systematic effect in our 
measurements that would produce a peak at this temperature. The same peak was reproduced at 
roughly the same height in three separate runs. We concluded that the peak is a real feature that 
probably coincides with the onset of premelting on the (0001) surface of ice at T « —12 C. This 
onset temperature agrees with that measured by p] , which is another premelting measurement using 
a surface preparation similar to that used in the current experiments. This temperature is also just 
below the large dip in ctq that we believe indicates a region of partial premelting (see Region II 
below). It is not obvious why the onset of premelting would result in a peak in ctq, however. For 
this we can only speculate some unknown many-body phenomenon at the surface. 

Our measurements suggest somewhat counterintuitively that ice growth at the nearby tempera- 
tures of —12 C and — 13 C could be qualitatively different under some circumstances, which could 
be investigated in other experiments. Another interesting note is that this is the first time the onset 
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of premelting has appeared in any measurement as a sharp feature. This suggests that additional 
exploration of the detailed growth dynamics of the (0001) surface of ice near — 12 C could be fruitful. 

Region II: [— 6 > T > — 11 C]. In this region we believe that premelting is partially developed on 
the (0001) surface. Our intended meaning here of "partially developed" is that premelting produces 
a partially nonfaceted surface that enhances nucleation and thereby results in the lower values of ctq 
measured, as was previously suggested by Kuroda and Lacmann Here, as in [T], we cannot define 
a detailed molecular model, since the surface molecular dynamics must be quite complex at these 
temperatures. One's normal picture of a static solid surface with a small number of admolecules 
is certainly too simplistic here. Nevertheless, our data show that the growth is still well fit by a 
nucleation model with aintrinsic = exp(— (To/c)- We therefore put forth a simple picture in which a 
partially developed premelted layer produces a dip in ag in this temperature region. This picture is 
supported by our measurements even if we cannot yet specify a molecular model of the ice surface 
dynamics. 

Region I: [—2 > T > — 5 C]. At these temperatures we suggest that the premelted layer has 
become fully developed, leading to an increase in ctq compared to that in Region II, again following 
[T]. Here our picture is that of a solid ice surface covered by a thick quasi- liquid layer (QLL), with 2D 
nucleation occurring at the ice/QLL interface. Because the ice/QLL interface is a smooth, faceted 
surface, cto is higher than that found on a partially premelted surface. To our knowledge, this is the 
first time that 2D nucleation at a solid/quasiliquid interface has been definitively observed. 

Concluding, we see that ice crystal growth data such as these can be used both as a probe of the 
temperature dependence of premelting and as a measure of the effects of premelting on ice crystal 
growth dynamics. We find it quite remarkable that the growth dynamics can be summed up so 
concisely by a single function (Tq (T) . This is true even going through the transition from no premelting 
to partially developed premelting to a fully developed quasiliquid layer, as described above. Although 
the equilibrium structure of the ice surface changes dramatically over this temperature range, as does 
the equilibrium vapor pressure, the functional form a = exp(— tro/f) remains unchanged as the ice 
growth is everywhere described by simple 2D nucleation; essentially all that changes is the critical 
supersaturation parameter ao{T). 

This work was supported in part by the California Institute of Technology and the Caltech- 
Cambridge Exchange (CamSURF) program. 
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